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Abstract 

I show that a usual propagator cannot be defined for the pseudo-diffusion equation of the Q- 
functions. Instead, a forward-backward propagator is defined, which motivated a generalization 
of Cahill-Glauber interpolating operator. An algorithm is also given for squeezing Q functions 
directly, using one-dimensional diffusion propagators. 



! 1 Introduction 

<n : 

In previous papers [0, § , it was shown that the Q functions for an operator A (see definition in Eq. 
(|l3|) below) obeys the following partial differential equation 



V(p,q;a) Q(A : p,q;a) 




Q(A : p, q; a) = , where a := e 



2y 



where y is the squeezing parameter, as defined in dl4|). (Note that a here is the inverse of the 

OS 



A := e~ 2y in jlj, ^], so that the roles of p and q get exchanged.) This equation describes how the 
Q functions Q(p,q;cr) get changed in phase space (p, q) as the squeezing parameter a is increased. 
If the Q function belongs to a density operator p, then Q(p : p, q; a) is a probability distribution 
I function, which remains normalized to 1 for all a, i.e. if its integral is carried out over the whole 

^ ■ phase space. Eq. (||) was called pseudo-diffusion equation Jl|, because (a) it resembles the 

diffusion equation in 2 dimensions Q, where the parameter a plays the role of time, and (b) 
the coefficients of and -^p in (Q) have opposite signs. (Since a is a monotonically increasing 
function of y, I shall use the time analogy when referring to either a or y, and hope that this will 
not lead to any confusion). Therefore, this equation describes a diffusive process in the q variable 
and an infusive one in the p variable for all a. In this way a thin distribution along the p-axis get 
continuously deformed into a thin distribution along the g-axis, as a is increased from to 00. 

The (symplectic) Fourier transform of the above partial differential equation ([l]), with respect 
to the phase-space variables (p,q), yields an ordinary differential equation in a: 



da + A\ k a 2 * 



Q(k,x;a) = , 



where the symplectic Fourier transform is defined by 

Q(k, x- a):= If ^e-^P-H Q(p, q - a ) . (3) 

It is easy to check that the " kernel " 

K{k, x- a, fi) := e-l^-^e-^-V- 1 -/^ 1 ) ; (4) 
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is a solution of Hence, the inverse Fourier transform of the products KQ yields the solution of 
the pseudo-diffusion equation ([[]), for any given initial Q function Q(p, q; p) with fi < a: 

Q(p,q;a) = jj^S*^K(k,x-a,ii)Q(k,x^). 



dkdx ^i\ X p-kg] 
2vr 



-K(k, x- o,n)JJ d ^L e -A*v'-U\ Q(p / j q >. M) . (5) 

Eq. @ shows that in order to get the solutions Q(p, q; a) we must perform two double integrations, 
first over (p',q') and then over (k,x). Therefore, it would have been much more efficient, if we are 
allowed to exchange the order of two double integrations, and carry out the integration over the 
variables (k, x) first: the resulting integral, if it existed, would have been equal to the Fourier 
transform of the kernel K: 

G,{p -p',q- q'; a, M ) := J J ^ ^V)-*^)] K ^ x . ^ ^ 

dkdx ,t[x(p-p')-fc( g -g')] i[ fc 2 ((T _ At)+a .2 ((T -i_ M -i )] 



2tt 



and would have enabled us to calculate the squeezed distributions (from unsqueezed or less squeezed 
ones) directly by a single double integration, as follows 

Q(p,q;o-)=JJdp'dq'G (p-p,q-q';o-,n)Q(p',q';n), for a > (j, . (7) 

Unfortunately, the integral (||) does not exist, as we shall see in section 3. Nevertheless, I shall give 
an algorithm for calculating squeezed Q functions which is based on seperation of variables and 
uses 1-dimensional diffusion propagators. 



Next, I show that it is more natural to define 2-sided propagators (f6) for the pseudo-diffusion 
equation ([!]); these propagators depend on three squeezing parameters, the present a, the past \i <a 
and the future A > a. Therefore, I shall refer to these 2-sided propagators as forward-backward 
(or future-past) squeezing propagators. A special case of these propagators connects Q functions 
to their Wigner counterparts. This makes it possible to give a new interpretation of the Wigner 
function. In turn, this interpretation makes it sensible to define formally generalized Q functions 
and Q operators which depend on two squeezing parameters, a p and a q . On a subdomain of these 
parameters, which can be expressed in terms of an ordinary squeezing parameter a and a "thermal 
parameter" r, the above generalized Q operators can be identified with thermalized and squeezed 
Wigner operators (I hope to give details elsewhere). Partial differentiation of these operators with 
respect to a and r yields the generalized pseudo-diffusion equation ([II]) and the diffusion equation 
(47), respectively. 

In this paper, I review in Sec. 2 basic properties of the Wigner and the Q functions. Then 
I explain in Sec. 3 why the Fourier transform of the kernel K cannot exist. In Sec. 4 I give an 
algorithm for calculating squeezed Q functions. In Sec. 5 I define the two-sided propagators for the 
pseudo-diffusion equation. In Sec. 6 I give the new interpretation of the Wigner function. In Sec. 7 

1 define and study generalized Q functions and Q operators. Finally, in Sec. 8 I give a summary. 

2 Wigner and Q functions 

The Wigner representation of an operator A is defined formally by Q 

/oo 
(q-x\A\q + x) e 2ixp dx = 2 Tr (AW(p, q) ) , (8) 
-co 
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where I use the round kets to denote eigenstates of the position operator, Q | x) = x \ x), to 
distinguish them from the number states \n). The operators (one for each (p,q)) 

/oo 
\q + x)(q-x\ e 2ixp dx , (9) 

were called Wigner operators by Dahl ||. Obviously, 

POO 

W(0,0) = / \x)(-x\ dx , (10) 



is the parity operator. Therefore, the Wigner operators can be interpreted as displaced parity 
operators [5-8], because 

W(p,q) =D(p,g)W(0,0)Dt(p, g ) = D(2p, 2g)W(0, 0) , (11) 

where 

D(p, q) := exp [i(pQ - qP)] = exp [-pq] exp [-iqP] exp [ipQ] (12) 
is called the displacement operator, because D(p, q) 

| ^ = e ip(a+«/2) | s + g ). Since W(0,0) is 
unitary and Hermitian, it follows immediately from ( |i~l| ) that also the displaced parity operators 
are unitary and Hermitian: W _1 (p, q) = WT(p,g) = W(p, q), so that W 2 (p, q) = I. This means 
that the W(p, g) are observables with eigenvalues ±1, as was pointed out by Dahl || and Bishop 
and Vourdas 0. 

The Wigner representation (|8|) yields functions of two variables, p and g, which may be looked 
upon as phase-space variables. These 'Wigner functions' have interesting properties and are useful 
for various calculations 0]. They are often referred to as quasi-probability functions, because they 
can take negative values even when A is a positive operator, A > 0, such as a density operator p. 

In contrast, the Q representation yield nonnegative functions for positive operators A: These 
functions are defined as follows [|], ||] 

Q(A:p,q;C) = (pq;C\A\pq;0=Tr(AU(pq;C)) , where U(p,q;() :=\ pq]Q(pq;(\ (13) 

are projection operators on the squeezed states \pq; £), which are defined by [||, |J 

\pq;C) = D(p,g)S(C)|0> , where C := ye iv (-00 < y < oo) (14) 

and |0) is the ground state of a specific harmonic oscillator, a|0) = 0. (i.e. a is the annihilation 
operator with a definite frequency ujq; In this paper, we set % = m = ujq = 1.) In © D(p,g) is 
the displacement operator (|l2|), which generates coherent states when applied to |0), and 

i(Ca^-CV)], (.==9^) (!5) 



S(C) = exp 



is the squeezing operator, where the squeeze parameter y vanishes in the coherent-state limit. 

If A is a density matrix p, then its Q function Q{p : p, q; £) can naturally be interpreted as 
a probability distribution. To emphasize this fact, the Q functions were denoted by P in Jl]. ||], 
instead of Q here. 

For simplicity, I shall from now on discuss only squeezings which are pure boosts, without 
rotation, i.e. with ip = in (|i~4|), and use the squeezing parameter a := e 2y instead of y. 
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3 Non-existence of a 1-sided Propagator 



I shall now explain why Go, the Fourier transform of the kernel K, can not exist. First, we note 
that the Go in (^), if it exists, would be a special case, for [i = A (!), of the following product of 
two Gi factors: 



G(p -p',q- q'; a' 1 - A" 1 , a - fx) := Gi(p - p'; a" 1 - A -1 ) G^q -q'-a-fj,), (16) 



where 



/c 
-I 



27r' 



Mp-P^e-^^- 1 -^ 1 ) 



Gi(q -q';cr- //) 



^(a- 1 - A" 1 ) 



•oo 27T 
1 



=e 



Note that Gi(x — x'; r) is the propagator of the one-dimensional diffusion equation 



only for A > a , (17) 



only for /i < a . (18) 



d_ 1 d 2 

dr 4 dx 2 



G-l(x-x';t) = 0, for t>0. 



(19) 



The integral (|l7|) exists only for cr < A, and then it yields the p-propagator, but it diverges for a > A. 
In contrast, the corresponding integral (18) exists only for a > fi, but diverges for a < fj,. Thus, 
their product fllh)) cannot exist simultaneously, if /j, = A. Consequently, the desired propagator Go, 
as a Fourier transform of K, cannot exist. 



4 An algorithm for squeezing Q functions 



To overcome the above difficulty we recall 0] that the solutions of the pseudo-diffusion equation 
(H) can be obtained by the method of separation of variables: Writing the solution as a product of 
two functions, Q(p,q;a) = 6(p,a)ip(q,a), where 6 depends only on p and a, and ip depends only 
on q and a, we get 








6(p-a) 



(20) 



Since the first term in ( po| ) depends only on q and a, while the second term in (20) depends only 
on p and a, we conclude that each of them must be equal to a function of a only, which we denote 
by /(cr). But I shall consider here only the case f(a) = 0, since it turns out that f(a) ^ yields 
no new solutions of ([!]). For f(cr) = 0, Eq. ( |20| ) yields the following two equations: 




4cr 2 dp 2 



e(p;a) 



' d 


i d 2 N 




~ 4 dq 2 f 


d 


i d 2 


a- 1 


4 dp 2 



771 " 7773 V>fec0 



0(p;<t) 



(21) 



(22) 
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where we used 



J^p = —-ji ga-i in (HI)- We see that t/j obeys a 1-dimensional diffusion equation in 



q, where a plays the role of time. Similarly, 9 obeys a diffusion equation in p, but with <r _1 playing 
the role of time. 

Let us for definiteness discuss the case a > 1 (the a < 1 case can be dealt with by similar 
arguments): For a > 1 we can obtain the solution ip(q, a) by using the propagator (18). The action 
of the propagator for a > 1 can also be written formally as a differential operator Yi(<7 — l,q), 
acting on the unsqueezed function ip(q, 1): 



ip(q, a) 



dq' 



1 



exp 



(a - 1) 8 2 



dq 2 



=: Yi(<7-1, • 



(23) 
(24) 



The propagator in ( |23| ) will blow up for a < 1, and the integral in ( |23| ) will not exist, even if ip(q, 1) 
is a strongly decaying function, such as tp(q, 1) = exp[— q 2 /a], for I — a < a < 1: This is because, the 
coefficient of q' 2 in the exponent in the integrand of (23) would be positive: — — = > 0. 

In contrast, if we /irsi evaluate the integral pffi ) using a > 1, i/ien we can continue the resulting 
function ifi(q,o~) to 1 — a < a < 1, as we shall illustrate below. 

This suggests the following algorithm for calculating the squeezed p factor: First, calculate the 
solution 9(p, a) of the diffusion equation ([22]) by using the p-propagator for 1/a > 1. Then continue 
the resulting solution to 1/a < 1. One could also calculate the solution a) directly for 1/a < 1, 
by using the power series expansion of the differential operator Yi(<r _1 — l,p): 



9(p,a~ l ) =Y 1 (a- 1 -l,p)9(p,l) = exp 



- 1: 



cP_ 
dp 2 



9(p, 1) . 



(25) 



As an example, let us squeeze the Q function of thermal light [p|, whose density matrix is given by: 



1 



Pth 



— —Y 

n+l n 



n) 



n=0 



n xn 



n+l 



(n\ , where n 



1 



30-1 



(26) 



where = hLo/ksT and n is the mean number of photons. The (coherent) Q function of (26) is 
factorizable: 



QiPth ■ P,q; 1) 



1 



n+l 



exp 



p 2 + q 2 
2(n + l) 



where each factor turns out to be proportional to the diffusion propagator: 



#7,1)=%,1) 



1 



y/n + 1 



exp 



2irGi(q,2n + 2) . 



(27) 



(28) 



Since the unsqueezed factors are the same, ip(q, 1 
by using the propagator for a > 1: 



2(n + l) 

= 9(q, 1), we only need to evaluate one integral 



ip(q,a) 



dq - 



1 



2vr Gx(q,a + 2n + l) 



cr + 2n + 1 



2vr / dq'G x (q-q',a-l)Gi (q', 2n + 2) 

2 

'! 

exp 



a + 2H+1 



(29) 



We see that i[)(q; a) is well defined for — (2n + 1) < a < 1. Since 9(p, 1) = ifi(p, 1), we get 9(p, a) = 
ijj{p,a~ l ) : so that the squeezed Q function becomes Q(pth '■ P,Q]&) = , 4 ) {p^o'~ l )'4){q,a). Thus, 
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the squeezing of Q{pth ■ Vi 1) required only a single integration, which in this specific case was 
carried out by noting that the successive action of two propagators is equivalent to the action of one 
propagator, whose 'time variable' is the sum of the individual times: (a — 1) + (2ra + 2) = a + 2n + 1. 



5 Forward- Backward Squeezing Propagators 



The products in fll6|) will exist and be solutions of the pseudo-diffusion equation, if we relax the 
condition fj, = A and demand instead only p < a < A. Actually, these products yield a different 
solution for each 4-tupel (p', q', A, /x). Since the heart operator W is linear, any superposition of 
these solutions is also a solution. In particular, if we fix the squeezing parameters \x and A and 
integrate only over p' and q', we get solutions of the form 



f(p,q;a 1 ,<x) = j j dp'dq' G(p - p' ,q - q'; A 1 - a 1 ,a - fj,) f(p' , q'; A for fi < a < A , 

(30) 

for any given function f(p,q;X , /x), provided that the integrals (30) exist. We see that G in 
(|30| ) is acting as a propagator, which provides a solution of the pseudo-diffusion equation (||) for 
the squeezing parameter a in the range /i < a < A. I shall therefore call these G functions two- 
sided or forward-backward squeezing propagators of the pseudo-diffusion equation (||), since the two 
squeezing parameters, p and A, lie on opposite sides of a. These G solutions have the proper limit, 
which one expects from a propagator, if a is approached from opposite directions: 

lim G{p-p ,q- q\\~ x -a~ l ,a-p) = 5{p - p') 5{q-q) . (31) 

An extreme case of the squeezing propagators (|l^) is obtained by choosing fi = and A = oo. These 
squeezing parameters correspond to the values — oo and +oo of the y = \ ln<7 variable, respectively: 

G(p — p' , q — q'; <7 _1 , a) = — exp —o'(p—p) 2 —a~ 1 (q — q') 2 , for < a < oo . (32) 

7T L J 

For the choice fi = and A = oo in (p0|), a can now take any positive value. Moreover, the square- 



root factors in the two propagators cancel out. Thus, with the propagator (32), the relation (|3 
becomes simply 

f(p,q-a-\a)= JJ^^eM-^P-p') 2 -cr- l (q-q') 2 }f(p',Q'-,0,0), for a > . (33) 



6 New Interpretation of the Wigner Function 

The Wigner representation of the projection operator H(p, q; a) is given by 

U(p, q-a) = ^Tr[n(p, q; a)W(p>, q')] W(p>, q') 

= 2ll^ex P [-o-(p-pr-o--\q-q') 2 ]W(p>,q>). (34) 

where we used 

TT[n(p,q;a)W(p',q')} = (p, q; a\W(p', q')\p, q; a) = (p,q;a\D(p' ,q')W(0,0)D^p' ,q')\p,q;a) 

= (p-p',q- q';o-\p -p,q' - q-a) = exp[-a(p - p') 2 - a~ l (q - q') 2 } . (35) 
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By multiplying by an operator A and taking the trace, we get a similar relation between 
the Q function and its Wigner counterpart [||: 

Q(A:p,q;a) = JJ ^ exp[-a(p - p') 2 - a~\q - q'f] W(A : p', q'), (36) 

Comparing this relation with (|33"1), we realize that the Q function is a propagated Wigner function 
via by the above special 2-sided squeezing propagator (|32j). This led me to define an interpolating 
function Q(A : p,q;a p l ,a q ) in Eq. ( |38[ ) below, which for a q = a p = a yields the usual squeezed 
Q function Q(A : p,q;a), whereas for a p = oo and a q = 0, it yields the corresponding Wigner 
function. Hence, we can interpret the Wigner function as a generalized Q function, which is 
squeezed to 'infinite future' y p = +00 (o~ p = e 2yp = 00) in the p variable, but anti-squeezed to 
'infinite past' y q = —00 (a q = e 2yq = 0) in the q variable. 



7 Generalized Q functions 



The above interpretation suggests the following definition of generalized Q functions which depend 
on two different squeezing parameters, a q and a p : Since the Wigner functions correspond to infinite 
past and infinite future, we can use the 2-sided propagators to bring them into finite past a q and 
finite future a p (or "two different points of the present"), as follows 



1 



Q{A:p,q;—,a q ) : 



Or, 



dp'dq' 



7T 



Gi(p - p'; a' 1 ) Gi (q - q'; a q ) Q(A : p', q'; 0, 0) 



(37) 



dp'dq' 



IT 



exp[-a p (p - p'f - a q \q - q'f] W(A : p',q') , (38) 



where W(A : p,q) = Q(A : p, q;0, 0) was used. Note that, for now, we assume a q ,a p > in (|38|), 
but a p may be larger or smaller than a q . Later, we shall extend the range of o~ q and a p . 

Instead of studying Q functions, it is more useful to define and study generalized Q operators, 
which yield Q functions by taking the trace Tr(^4Q). Similar to (|38|), we generalize the operator 
relation (34), as follows: 

Q(p,^;c r p 1 ,c r g) 



2 exp 
2 

f f dkdx 
J J ~2T 



1 / 1 d 2 d 2 ^ 

+ (Tq dq- 2 



a p dp 2 



W(p,q) 



(39) 



dp'dq' 



IT 



exp[-a p (p -p') 2 - a-\q - q') 2 } W(p',q') (40) 



■ exp 



• — ( — x 2 + o~ q k 2 ) + i(px — qk) 
4 cr p 



D(fc,x) 



(41) 



where (^) follows from ( p0| ) by substituting the Fourier transform of the Wigner operator ||: 



W(p,g) = - 



jj dM^ el[px _ qk]n{k ^ x) 



(42) 



So far we assumed that a p and a q are positive. However, we can formally extend the definition 
(|39|) to real (and even complex) values of a p and a q . 

From (^) or ( ]4l| ) we see immediately that the Q operators (and hence also the generalized Q 
functions via the trace) satisfy the following partial differential equation 







1 / da q d 2 



1 da p d 2 



dfj, 4 I dfj, dq 2 a 2 d[i dp 2 



o-p(n) 



(43) 
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if o q and o p are varied along a curve (o q (fj,),o p ([i)) in the (o q ,o p ) plane. We see that 



yields 



a generalized pseudo-diffusion (diffusion) equation, if and have the same (opposite) signs. 
In particular, on the subdomain of M 2 , defined by 



-i 



OT 



where o > , r > — 1 . 



j p — o~t " , o~ q = o~t , wnere o > u , r > — i . (44) 

the operator Q(p, q; to~ 1 , to) yields a generalization of the Cahill-Glauber interpolation operator 
[10 1, which corresponds to zero squeezing (cr = 1): 



Q(p,q;TO 1 ,rcr) 



dkdx 
~27 



exp 



(<t x + crfc ) + i(px — (/&) 



D(fc,x) . 



(45) 



We see that the parameter s of [10] corresponds to our — r here. I shall call r the thermal variable, 
since Q(0, 0; r, r) is equal to the thermal density operator p t h in ([26|), if we make the identification 



r = coth(/3/2). Hence, Q(0,0;r, r) is equal to the thermalized parity operator of [11] and to the 
generalized parity operator of which was introduced in order to obtain the parity operator 
W(0, 0) as a limit of trace-class operators. 

From equation (|43|) we see that Q(p, q; rcr -1 , to) obeys the following two interesting partial 
differential equations: 



^(p,<?;T,cr) Q(p,q;To 1 ,to) 



d 


T 


/ d 2 


i d 2 


do 


~ 4 


\dq 2 ' 


o 2 dp 2 


d 




1 o ^ 


i d 2 

o dp 2 


8t 


-il 


V dq 2 



Q(p,q-TO- l ,To) = 0, (46) 
Q(p,q-TO- l ,To) = 0, (47) 



if o is varied along 'isothermal curves' (o q /o p = t 2 = const.), or r varied along the 'isosqueeze 
curves' (o q o p = o 2 = const.), respectively. 

Eq.(|46|) is a generalized pseudo-diffusion equation, since ^(p, q; r, o) now depends also on r and 
reduces to the original heart operator in (Q) for r = 1. For r = equation (^6|) tells us that the Q 
operator is not affected by squeezing: Q(p, q;0, 0) = 2W(p, q) for all o. 

Eq.(47) is a diffusion equation in 2 dimensions Q, with different diffusion constants in the p 
and q directions, if o ^ 1. 



8 Summary 

I reviewed the 3-steps Fourier-transform procedure for calculating squeezed Q functions Q(A : 
p, q; o) from given unsqueezed or less squeezed ones Q(A : p, q; fi), where [i < o 0: First calculate 
Q(A : k,x;fi), the Fourier transform of Q(A : p,q;fi), then multiply Q by the special kernel K 
and finally take the inverse Fourier transformation of the product KQ and you get Q{A : p,q;o). 
Instead, I gave an algorithm for squeezing Q functions directly, by using one-dimensional diffusion 
operators. 

I also explained why one-sided propagators cannot be defined for the pseudo-diffusion equation 
(H) and showed that two-sided squeezing propagators (|l~6|) are more appropriate for this equation. 

I noted that the Q functions are related to their Wigner counterparts by the extreme two-sided 
propagator (^), and concluded that the Wigner functions can be looked upon as generalized Q 
functions Q(A : p, q; o~ , o q ), which are squeezed forwards (o p = oo) in p variable and backwards 
(o q = 0) in the q variable. 

These Q operators were defined formally in ( |39| ) in terms of Wigner operators. On a subdomain 
of the parameters (o p ,o q ), which is defined in terms of r and o in (|J|), these Q operators yield 



S 



the 'thermalized and squeezed Wigner operator' Q(p, q; t<t _1 , to-), which yields the interpolating 
operator of Cahill and Glauber ju]] for zero squeezing (a = 1). However, we emphasize that even 
Q(p, q; to" -1 , Ta) for a > is less general than our generalized Q operator Q(p, q; o~ , a g ): For 
example, Q(p, q;0, a q ) does not correspond to a thermalized and squeezed Wigner operator. 

I also showed that Q(p, q; to~ x , Ta) obeys the generalized pseudo-diffusion equation ([46]) for 
constant r, and the 2-dimensional diffusion equation ( f47| ) for constant a. This diffusion equation 
(47) explains intuitively why each of the P, the Wigner and the Q distributions, become smoother 
than its predecessor, as r is increased from r = —1 to r = to r = 1, for any a = const. 

This smoothing process is expected to continue as r is increased from 1 to oo. Indeed, it can 
be shown that the generalized Q functions Q(p : p, q; TO~ l ,Ta) = Tr [pQ(p, q; t<t -1 , rcr)] of density 
matrices p are nonnegative for r > 1, and thus they yield probability distributions and not merely 
quasi-probability distributions. Consequently, their Wehrl entropy [|l^] 

S(T,a)\=-jj^—Q(p-.p,q]Ta' l ,Ta)\nQ(p:p,q]T(T' 1 ,Ta), for r>0, (48) 

is well defined and should increase if we vary both r and a simultaneously in such a way that 
both o q and a p increase. Hence, further clarification of the physical meaning of the observable 
Q{p, q\ro~ x ,t<t) would be useful. 
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